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ABSTRACT 



We present X-ray observations of the ultraluminous infrared galaxy, UGC 
5101, thought to contain a buried active galactic nucleus (AGN) based on ob- 
servations in other wavebands. We detected an absorbed hard component at >3 
>■ ■ keV, as well as soft emission in the energy range 0.5-2 keV. The soft X-ray com- 

ponent, possibly due to a modestly dust-obscured, extended starburst, has an 
' absorption-corrected 0.5-2 keV X-ray luminosity of Lx (0.5-2 keV) = 1.2 x 10^^ 

^ ■ ergs s~^. The 0.5-2 keV X-ray to infrared luminosity ratio is a factor of ~5 lower 

than typical values for a normal starburst, suggesting that this extended star- 
burst is unlikely to be energetically dominant in UGC 5101. The most plausible 



o 

Q^' origin of the absorbed hard component is the putative buried AGN. The 6.4 keV 

O ■ Fe Ka emission line has a modest equivalent width (~400 eV), suggesting that 

' this hard component is direct emission from the AGN, rather than a scattered 

ci • component. The absorption-corrected 2-10 keV X-ray luminosity of the buried 

AGN was estimated to be Lx(2-10 keV) ~ 5 x 10"^^ ergs s~^ The intrinsic 2-10 
^ ■ keV X-ray luminosity and the 2-10 keV X-ray to infrared luminosity ratio are 

■ both comparable to values measured for Mrk 463, a Seyfert-2 galaxy of similar 

infrared luminosity. 

Subject headings: galaxies: nuclei — X-ray: galaxies - galaxies: individual (UGC 
5101) 
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1. Introduction 

Ultraluminous infrared galaxies (ULIRGs) detected by the IRAS all-sky survey radiate 
large, quasar-like luminosities (> a few x 10'^^ ergs s~^) as infrared dust emission (Sanders et 
al. 1988); therefore they must consist of extremely powerful energy sources that arc obscured 
by dust. ULIRGs have also been found at large distances {z > 0.3), where they have been 
used to trace the dust-obscured star formation rate, dust content, and metallicity of the early 
universe (Blain et al. 1999). However, it remains unclear whether ULIRGs are powered by 
dust-obscured active galactic nuclei (AGNs) or starbursts. 

It is now widely accepted that the energetically dominant components in the majority of 
ULIRGs are compact (<a few 100 pc), highly dust-obscured nuclear energy sources, rather 
than extended (~kpc), weakly obscured starburst activity (Soifer et al. 2000). It is essential 

to determine whether these compact energy sources consist of AGNs and/or compact star- 
bursts. Approximately 30% of ULIRGs classified optically as Seyferts are known to possess 
energetically important AGNs surrounded by dust with the geometry of a torus (Veilleux, 
Kim, & Sanders 1999a; Veilleux, Sanders, & Kim 1999b). However, the AGNs expected to 
reside in the majority of non-Seyfert ULIRGs may be deeply buried in dust along all lines 
of sight. It is necessary to determine the energetic importance of such buried AGNs in order 
to understand the true nature of ULIRGs. 

To estimate the relative energetic contributions of an AGN and a starburst, it is nec- 
essary to utilize signatures that are produced by only one of these processes. One such 
signature is the polycyclic aromatic hydrocarbon (PAH) emission found in the infrared 3-15 
lim wavelength range (AUamandola et al. 1989). PAH emission is usually observed in a 
normal starburst, but not detected in a pure AGN (Moorwood 1986; Roche et al. 1991). 
Thus, PAH emission can be used to detect starburst emission, by disentangling from AGN 
emission, in an AGN/starburst composite, such as a ULIRG, and PAH luminosity traces the 
starbutst's absolute magnitude (Mouri et al. 1990; Genzel & Cesarsky 2000; Imanishi 2002). 
The 3.3 fim PAH emission feature and dust absorption features found in infrared 3-4 fim 
spectra have been used to detect candidates for non-Seyfert ULIRGs possessing energeti- 
cally important buried AGNs (Imanishi & Dudley 2000; Imanishi, Dudley, & Maloney 2001; 
Imanishi & Maloney 2003). 

One example of this type of object is UGC 5101 (z = 0.040). It has an infrared (8-1000 
Hm) luminosity of Ljr ~3 x 10"^^ ergs s~\ a cool far-infrared color {IRAS 25/um-to-60/xm 
flux ratio of <0.2), and is optically classified as a LINER, based on a systematic spectral 
classification of galaxies (Veilleux et al. 1995). Imanishi et al. (2001) found the observed 
infrared 3-4 yum spectrum of UGG 5101 to be best explained by the presence of a powerful 
buried AGN. Imanishi & Maloney (2003) argued that the nuclear energy source in UGG 5101 
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is more centrally concentrated than nuclear gas and dust, as is expected for a buried AGN, 

rather than a normal starburst, which shows a mixed dust geometry (Soifer et al. 2000). 
Radio VLBI observations (Smith et al. 1998; Lonsdale et al. 2003) and a near-infrared high 
spatial resolution image obtained with HST (Scoville et al. 2000) also detected compact, high 
surface-brightness emission, which was interpreted as an AGN by these authors. Therefore, 
there is mounting evidence for the presence of a buried AGN in this galaxy. The detection 
of strong, absorbed 2-10 keV X-ray emission can provide compelling evidence for a powerful 
buried AGN, because an AGN is a much stronger 2-10 keV X-ray emitter than a starburst. 
Previous X-ray observations with Chandra failed to find clear evidence for X-ray emission 
from the putative AGN (Ptak et al. 2003). Here, we report analysis of high quality X-ray 
data of UGC 5101 obtained with X MM and Chandra. The high spectral sensitivity of XMM 
at 2-10 keV was fully utilized to detect absorbed hard X-ray emission from the putative 
buried AGN. Chandra data arc used to investigate the X-ray morphology, benefiting from 
its high image quality, as well as to check the overall spectral consistency with the XMM 
spectrum. Throughout this Letter, we use a cosmology of Hq = 75 km s~^ Mpc^^, flu = 
0.3, and Qa = 0.7. 



2. Observations and Data Analysis 

An XMM observation of UGC 5101 was performed with the EPIC cameras (PN, MOSl, 
and M0S2) in full-frame mode on 2001 November 12. We used standard data analysis pro- 
cedures with the XMM Science Analysis Software (SAS version 5.4), together with standard 
software packages (FTOOLS 5.2 and XSPEG 11.2). Time intervals with high background 
count rates were removed from the events files, and we included only events corresponding 
to patterns 0-4 for PN and 0-12 for MOS, as suggested in the SAS Handbook. The resulting 
exposure times, after screening for time intervals with high background, were 25.1 and 33.1 
ksec for PN and MOS, respectively. We used the most recent calibration files retrieved from 
the XMM web page and response matrices created with SAS for our data analysis. 

We extracted a source spectrum using a 22.5 arcsec radius circular region. A back- 
ground spectrum was extracted from a nearby source-free circular region with a radius of 
60 arcsec. We binned the spectra so that each energy channel contained a minimum of 20 
counts. An absorbed hard X-ray component was clearly detected only in the PN spectrum, 
which contains the most counts. Therefore, we used the PN spectrum to obtain rough mea- 
surements of important parameters. However, to improve the statistics, final spectral fitting 
was performed using the combined PN, MOSl and M0S2 spectra. We allowed the relative 
normalizations of the cameras to vary in our spectral fits as there is some uncertainty in the 
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relative flux calibration of the three cameras. 

A Chandra observation was performed on 2001 May 28 with the ACIS-S3 back-illuminated 
CCD chip. The data were reprocessed using CIAO 2.2.1 and CALDB 2.17. An effective ex- 
posure time of 47.8 ksec was obtained after discarding high background intervals. 

A Chandra spectrum of the nucleus was extracted from a circular region with a radius 
of 8 arcsec. A background spectrum was extracted from an annular region around the 
nucleus and subtracted from the nuclear spectrum. The spectrum was binned to facilitate 
the detection of emission-line features. The C-statistic was used to fit the spectrum since 
the number of events in each bin is small. 



3. Results 

Figure 1 shows a 0.5-2 keV Chandra image. The source is extended beyond the PSF of 
Chandra (~1 arcsec diameter for >85% encircled energy). 

Figure 2a shows the XMM EPIC spectrum of UGC 5101. In addition to the soft emission 
dominating the 0.5-2 keV X-ray range, an absorbed hard component at >3 keV and emission 
line at ~6.5 keV were clearly seen. The hard component is likely to originate in the putative 
buried AGN, and so we fitted this component with a power-law model. The emission line 
at ~6.5 keV is hkely to be the neutral Fe Ka line at 6.4 keV, and was fitted with a narrow 
Gaussian profile. For the soft component, a 0.5-1 keV thermal emission model (mekal model 
in XSPEC) from a starburst was first attempted to fit the soft X-ray spectrum, particularly 
the apparent spectral peak at 0.7-1 kcV. However, significant residuals were present in this 
model fit. Therefore, we added (A) higher temperature thermal emission (5-10 keV) from 
hot gas and X-ray binaries in a starburst, or (B) a second power-law component, which could 
be due to Thomson scattering of the power- law emission from the AGN. Table 1 summarizes 
the best-fit parameters for these models. To estimate the power of the putative buried AGN, 
the most important measurement is the absorption-corrected 2-10 keV X-ray luminosity of 
the hard component, Lx(2-10 keV). Both models result in a similar Ljss:(2-10 keV) of ~5 x 
20^2 gj,gg g-i g^j^j have a 6.4 keV Fe Ka line of modest equivalent width (EW6.4 ~ 400 eV). 
We varied some of the model parameters within reasonable ranges (e.g., a power-law photon 
index of F = 1.0-2.0, a temperature between 5-9 keV for the high-temperature thermal 
emission, and abundances of 0.1-0.4 solar for the thermal emission); however, these changes 
result in values of Lx(2-10 keV), which are consistent within a factor of <2. 

Figure 2b shows the Chandra ACIS spectrum. The solid line overplotted for comparison 
is based on model A (in Table 1) from the XMM EPIC spectrum, with normalization of the 
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continuum and of the Gaussian component representing the Fe Ka line as free parameters 
and were determined from the Chandra ACIS spectrum. The overall spectral shape is very 
similar to the XMM spectrum. In the soft energy band, several emission lines are seen: strong 
He-like Si (~1.9 kcV) and He-likc S (~2.5 kcV). We detected the absorbed hard component, 
but the amount of absorption of this component is not well constrained from the Chandra 
spectrum alone, because of the paucity of photons at the highest energies (>7 keV). 

The Fe Ka emission feature at 6.4 keV is also detected in the Chandra spectrum. To 
estimate its equivalent width, we fitted the Chandra spectrum in the 4-7.5 keV band to 
minimize uncertainties in the modeling of the continuum. Model A (in Table 1) was applied 
to the continuum. The equivalent width was estimated to be EW6.4 = 230^12 eV (Ax^ = 
2.7). This value is similar to the estimate based on the XMM EPIC spectrum (Table 1). 
The extremely large equivalent width (5.9^3^2^ keV, A^^ = 4.605) reported by Ptak et al. 
(2003) was not confirmed with our analysis. 

4. Discussion 

4.1. Energetic Importance of Extended Starburst Activity 

UGC 5101 shows starburst activity with spatial extension of a few kpc (~4 arcsec) in 
the east- west direction (Gcnzcl ct al. 1998), with dust obscuration at a level similar to that 
of the nearby prototypical starburst galaxy M82 in both mixed-dust and foreground- screen 
dust models (McLcod ct al. 1993; Sturm et al. 2000). The extended morphology seen in the 
Chandra 0.5-2 kcV image (~4 arcsec in the cast-west direction; Fig.l) and the detection 
of emission lines in our spectra (Fig. 2) suggest that thermal emission from the extended 
starburst contributes significantly to the soft X-ray emission below 2 keV. Assuming that all 
of the soft X-ray emission is of thermal origin (model A in Table 1), the absorption-corrected 
0.5-2 keV soft X-ray luminosity is Lx(0.5-2 keV) ~ 1.2 x 10^^ ergs s~^, which results in 
Lx(0.5-2 keV)/LiR ~ 4 x 10"^ 

In a starburst galaxy, the 0.5-2 keV X-ray to far-infrared (40-500 /im) luminosity ratio 
is typically Lx(0.5-2 keV)/LFiR ~ 2 x 10"^ (RanaUi, Gomastri, & Setti 2003). As the far- 
infrared (40-500 /jLUi) and infrared (8-1000 //m; Sanders & Mirabel 1996) luminosities do 

not differ very much in an ordinary starburst, the Lx(0.5-2 kcV)/LiR ratio in UGC 5101 is a 
factor of ~5 lower than the value expected for a normal starburst. If some fraction of the soft 
X-ray emission originates in scattering of the power-law AGN component, this discrepancy is 
even larger. In UGC 5101, the ratio of starburst-tracing 6.2 fim PAH luminosity (measured 
using large-aperture spectroscopy with ISO; Spoon et al. 2002) to infrared luminosity, is 
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also lower by a factor of ~4 than that of a normal starburst (Imanishi & Maloney 2003). As 
the properties of the extended starburst in UGC 5101 are not expected to differ significantly 
from a normal starburst, these results both suggest that the modestly obscured extended 
starburst is unlikely to be the dominant energy source in UGC 5101. 

4.2. A Buried Energy Source 

As the hard X-ray component is strongly absorbed, it cannot originate in X-ray binaries 
from the modestly obscured extended starburst (Persic & Raphaeli 2002). The most plausible 
origin is an AGN deeply buried in gas and dust at the core. The modest equivalent width 
of the 6.4 keV Fe Ka line suggests that direct X-ray emission from the AGN is detected, 
because in physically plausible AGN geometries, a scattered X-ray component is expected 
to show a very large (>1 keV) equivalent width (Matt, Guainazzi, & Maiolino 2003). If 
the absorbed hard component is direct emission, we find an absorption-corrected 2-10 keV 
luminosity for the buried AGN of Lx(2-10 keV) ~ 5 x 10^^ ergs (§ 3). In principle, the 
modest EW6.4 of the detected hard component could be produced in a scattering situation, 
if it is Thomson scattered by highly ionized gas. If this is the case, the intrinsic Lx(2-10 
keV) is even higher than our estimate above. 

The Lx(2-10 keV)/LiR ratio for UGC 5101 is ~2 x 10"^ Lx(2-10 keV) values have 
been measured for some Seyfert-2 infrared luminous galaxies (Lir > a few xlO^^ ergs s~^) 
where the detected 2-10 keV emission is interpreted as a direct component, because of the 
modest equivalent widths of the 6.4 keV Fe Ka emission (<1 keV). Examples include IRAS 
20460+1925 (Ogasaka et al. 1997), IRAS 23060+0505 (Brandt et al. 1997), PKS 1345+12 
(Imanishi & Ueno 1999), and Mrk 463 (Ueno et al. 1996). All are thought to be powered 
by AGNs, with no detectable starburst activity in the infrared (Genzel et al. 1998; Veilleux, 
Sanders, & Kim 1999b; Imanishi & Dudley 2000; Imanishi 2002). The Lx(2-10 keV)/LiR 
ratios of these known AGN-powered, Seyfert-2, infrared luminous galaxies are between 0.002 
and 0.02, with Mrk 463 having the lowest value (Lx(2-10 keV) ~ 4 x 10^^ ergs s"^; Ueno et 
al. 1996, Lir ~ 2 x 10^^ ergs s~^). A smaller value of Lx(2-10 keV)/LiR is expected for a 
buried AGN than for a Seyfert-2 AGN surrounded with dust in a toroid geometry, because 
in the former case, a larger fraction of the photons emitted from the AGN are absorbed by 
dust and re-emitted in the infrared. Thus, the buried AGN in UGC 5101 is likely to be as 
powerful as that in Mrk 463, in terms of both absolute 2-10 keV X-ray luminosity and ratio 
to infrared luminosity. 

Regardless of whether the primary energy source of UGC 5101 is a buried AGN or a 
starburst, the large infrared luminosity is produced by UV emission. The contribution from 
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the detected X-ray-emitting buried AGN to the total energy output in the infrared depends 
highly on the intrinsic UV to 2-10 keV X-ray luminosity ratio; this is poorly understood and 
could vary substantially among the population of AGNs. The buried AGN could account 
for the bulk of the infrared luminosity in UGC 5101 if its intrinsic UV to 2-10 keV ratio is 
similar to the AGN in Mrk 463. 

By applying a template of infrared spectral energy distribution of a classical, moderately 
infrared luminous AGN and a normal starburst to UGC 5101, Farrah et al. (2003) argued 
that infrared dust emission flux from UGC 5101 at A < 10 fim is dominated by an AGN, 
while that at A = 10-1000 /xm is dominated by a starburst. This conclusion comes primarily 
from the known fact that a classical AGN generally shows a warmer infrared color than a 
normal starburst (de Grijp, Miley, & Lub 1987). However, ULIRGs are known to have a 
larger column density of obscuring dust than classical, less infrared luminous AGNs (Spoon 
et al. 2002; Imanishi & Maloney 2003), and thus infrared emission from the buried AGNs 
in ULIRGs should be cooler than that from classical AGNs. Furthermore, a significant 
fraction of dust grains in UGC 5101 was found to be ice-covered (Spoon et al. 2002; Imanishi 
& Maloney 2003), which also causes the infrared color to be cooler than bare dust grains 
(Dudley, Imanishi, & Maloney 2003). Therefore, a much larger fraction of the infrared 
dust emission from UGC 5101 at >10 can be produced by the buried AGN than that 
estimated by Farrah et al.. Indeed, if an observed excess in the 60 /im luminosity function at 
the highest luminosity end, over the extrapolation from lower luminosity, is caused by AGN 
contribution (Takeuchi et al. 2003), then it is suggested that AGN-powered dust emission 
can contribute significantly to the luminosities of ULIRGs at as long as 60 fj,m. These factors 
should be taken into account to correctly understand the origin of strong dust emission from 
UGC 5101 at A > 10 yum. 

Y.T. and N.A. are supported by the Japan Society for the Promotion of Science. The 
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Table 1. Fitting results of the XMM data of UGC 5101. 



Model 
(1) 


xVd.o.f. 
(2) 


kT 
[keV] 
(3) 


Nh 
[10^3 cm-2] 

(4) 


LTH(0.5-2keV) 
[10^0 ergs s-i] 
(5) 


LpL(2-10keV) 
[10^2 ergs s~^] 
(6) 


EW 
[eV] 
(7) 


A 
B 


45.7/46 
50.5/46 


0.7l°| 
0.6l°i 


14+^ 
121^ 


12 
3 


7 
6 


410li° 



Note. — Col. (1): Adopted best-fit models. Model A is absorbed power law (F = 1.8 fixed) + 0.7 
keV thermal + 10 keV (fixed) thermal + narrow Gaussian (6.4 keV Fe Ka). Model B is absorbed 

power law (F = 1.8 fixed) + 0.6 keV thermal + second power law (F = 1.8 fixed; scattered 
component) + narrow Gaussian. In both models, the Galactic absorption is included. For thermal 
components, a mekal model, a fixed abundance of 0.5 solar, and plasma density of 1 cm"^ were 

adopted. Col. (2): Reduced value. Col. (3): Temperature for soft thermal emission. The 
uncertainty is for 90% confidence level for one parameter of interest (Ax^ = 2.7) throughout this 
Table. Col. (4): Absorption for the main power law component. Col. (5): Absorption corrected 
0.5-2 keV luminosity for thermal emission. In model A, both 0.7 keV and 10 keV thermal 
components were added. Col. (6): Absorption corrected 2-10 keV luminosity for the primary 
power law component. Col. (7): Equivalent width of the 6.4 keV Fe Ka emission line, relative to 

the main power law component. 
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Fig. 1.— A Chandra image of UGC 5101 in the 0.5-2 keV X-ray band. 
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Fig. 2.— (a): An XMM EPIC (PN, MOSl, and M0S2) spectrum of UGC 5101, with model 
A in Table 1 overplotted. In the upper panel, the higher plots are a PN spectrum (which 
contains the most counts) , and the lower plots are MOS spectra. Residuals of the data from 
the best-fit model are displayed in the lower panel, (h): A Chandra ACIS spectrum of the 
nuclear region. The model A for XMM EPIC data is overplotted as a solid line. The lower 
panel shows the ratio between the data and the model. 



